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ABSTRACT 
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The  Hamiltonian  for  a  system  of  interacting  electrons  in  the  presence  of  a 
barrier  is  transformed  into  one  which  can  be  separated  into  three  parts »  Two  of 
the  parts  describe  the  electrons  on  the'^right^and -^eff'^sides  of  the  barrier, 
while  the  third  is  a  transition  term  which  allows  tunnelling  through  the  barrier. 
It  is  shown  that  it  is  not  reasonable  to  make  the  '*rightn^ and  '*ieft'1f^ hand 
Hamiltonians  commute,  but  the  effect  of  this  failure  is  easily  accounted  for  in 
perturbation  theory.  In  addition,  corrections  to  the  transition  operator  arising 
from  inter-electronic  interaction  are  written  down  and  discussed  briefly.  They 
may  be  both  quantitatively  and  qualitatively  important, 
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TUNNELLING  FROM  A  MANY  PARTICLE  POINT  OF  VIEW 


I. 

Experiments  in  which  electrons  are  observed  to  tunnel  through  insulating  layers 

(1) 

have  played  an  increasingly  important  role  in  recent  times.  The  tunnelling  pro¬ 
cess  has,  of  course,  been  well  understood  since  the  discovery  of  quantum  mechanics. 

(2) 

Bardeen  has  very  recently  pointed  out,  however,  that  in  the  case  of  a  many- electron 
system,  the  tunnelling  process  is  very  conveniently  regarded  as  in  time-dependent 
perturbation  theory.  In  other  words,  it  is  seen  as  a  transition  process  from  one  set 
of  nearly  Stationary  states  to  another. 

Let  us  keep  in  mind,  to  be  definite,  a  system  of  interacting  electrons  enclosed 
in  cubical  boxes  separated  by  a  high,  narrow  barrier.  It  will  be  a  convenient  simpli¬ 
fication  in  the  later  development  if  we  suppose  that  the  boxes  and  the  barrier  are 
symmetrically  disposed,  as  in  Fig.  1. 


Fig  I. 


Barrier  Potential  V(x)  vs.  x. 
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The  most  convenient  formulation  of  the  problem  will  involve  showing  that  the 
complete  Hamiltonian  of  the  many  electron  system,  including  the  insulating  barrier,  can 
be  split  into  the  form, 

H  =  Hp  +  Hl  +  T  (1) 

is  the  Hamiltonian  of  the  electrons  to  the  "right"  of  the  barrier,  H^  the 
Hamiltonian  for  the  left-hand  electrons,  and  T  is  the  operator  which  makes  possible 
the  transition  of  electrons  from  left  to  right.  It  is  perhaps  clear  intuitively 
from  the  work  of  Bardeen  that  such  a  formulation  of  the  problem  exists,  and  indeed, 

,  x  <3) 

several  papers  have  appeared  in  which  the  Hamiltonian  (1)  was  taken  for  granted. 

It  is  the  purpose  of  this  note  to  show  in  what  sense  (1)  is  actually  correct. 
Perhaps  this  is  not  an  entirely  empty  exercise  since  an  instructive  difficulty  has 
to  be  overcome.  Furthermore,  it  is  possible  to  discuss  in  detail  the  effect  of 
interactions  on  the  above  formulations,  which  is  always  difficult  to  do  in  an 


Here  H^ 


intuitive  fashion 
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The  dilemma  which  arises  is  as  follows s  In  defining  the  "right-hand"  and 
"left-hand"  Hamiltonians,  one  naturally  needs  to  define  right  and  left-hand  states. 

Now,  first  of  all,  these  wave-functions  have  to  be  in  the  domain  of  the  original 
Hamiltonian  operator,  H.  In  the  concrete  model  we  have  adopted  this  simply  means 
that  they  must  vanish  at  the  sides  of  the  boxes  jx)  ~  a,  as  well  as  being  twice 
differentiable.  Secondly,  we  would  like  the  left-hand  states  to  be  complete  on 
the  left-hand  side.  This  requirement  insists  that  should  an  electron  flow  from  a 
wire  into  the  box  to  the  left  of  -  b  ,  its  wave- function  could  be  expressed  in  terms 
of  left-hand  states  only.  Thirdly,  one  should  like,  if  possible,  that  all  left-hand 
states  be  orthogonal  to  all  right-hand  states,  in  order  that  should  commute  with  H^. 
Finally,  the  states  have  to  form  a  complete  set  when  taken  all  together. 

If  one  could  produce  such  a  set  of  (single-electron)  wave  functions  0j^(x) 
and  ,  it  would  be  a  simple  matter  to  expand  the  annihilation  operator  as 

(4*  >  =  £  &r}*>  +  4  Ktx>  (2) 

One  would  then  replace  the  creation  and  annihilation  operators  in  the  second- 
quantized  form  of  H, 

j~f  c  l/lx>  +  if 4^  w 

by  the  expression  (2) ,  H^  could  be  identified  as  that  part  of  H  containing  only 
"a"  operators,  H  would  contain  only  "b"  operators, and  T  would  be  the  remainder. 

Unfortunately,  no  such  states  appear  to  exist.  For,  should  all  of  the 
(p^( X )  vanish  identically  for  positive  values  of  the  coordinate,  then  either 

they  do  not  form  a  complete  set  for  negative  x,  or  they  are  too  singular  at  x  =  0. 

And  if  they  have  support  in  the  right-hand  region,  they  cannot  be  orthogonal  to 

the  (ft  (x)  should  these  functions  be  complete. 

* 
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We  are  in  consequence  forced  to  compromise  the  requirements  on  the  wave 
functions.  There  are  two  approaches  which  can  be  taken.  The  first  possibility 
is  to  maintain  the  requirement  of  mutual  orthogonality,  but  to  abandon  the  re¬ 
quirement  of  separate  "right"  and  "left"  completeness.  Instead,  the  right  and 
left  states  will  be  made  "as  complete  as  possible"  in  their  respective  regions. 

If  this  is  done,  a  mathematically  convenient  formalism  results,  since  the 
operator  T  can  be  treated  systematically  by  ordinary  perturbation  theory.  We 
shall  show  below  a  simple  method  for  achieving  this  result.  One  must  pay  a  price, 
however.  In  practice,  the  model  which  we  are  studying  is  not  really  of  interest 
by  itself.  We  also  need  a  method  of  hooking  up  wires  and  batteries  connecting  the 
right-hand  electrons  with  those  on  the  left.  Now  wires  and  batteries  are  cumber¬ 
some  to  treat  quantum  mechanically,  to  say  the  least.  It  is  essential  to  have 
a  simple  and  convincing  method  of  mathematically  treating  the  classical  part  of 
the  apparatus.  These  methods  usually  involve  introducing  electrons  into  the 
right-hand  states  and  taking  them  out  of  the  left-hand  ones.  However,  the  wires 
introduce  electrons  into  the  right-hand  side.  If  the  right-hand  states  are  not 
complete,  we  have  to  take  into  account  the  possibility  that  electrons  flowing  into 
the  right  side  partially  go  into  the  left-hand  states .  Methods  may  be  advanced  to 
get  around  this  difficulty,  but  they  seem  rather  contrived.  The  approach  does  con¬ 
tain  the  possibility  of  conveniently  studying  higher  order  effects,  which  is  be- 

W 

coming  fashionable. 

Another  approach  is  to  abandon  the  orthogonality  condition,  requiring  instead 
the  states  to  be  "as  orthogonal  as  possible",  and  to  keep  the  separate  completeness 
of  the  right  and  left  states.  In  that  case,  the  commutator  Or*  Hl]  is  an 
operator  of  the  order  of  magnitude  of  T, 
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in. 


In  the  symmetrical  model  which  we  are  considering,  a  particularly  simple 

set  of  states  springs  to  mind.  The  exact  eigenstates  of  the  single-particle 

Hamiltonian  may  be  called  (x) ,  Or(x) .  They  have  even  and  odd  parities  under 

the  reflection  of  the  x-coordinate,  respectively.  These  states  occur  in  pairs,  the 

energy  of  the  even  (odd)  state  being  £  (T)T  ,  to  an  excellent  approximation, 

k  Kk 

where  T^  happens  to  be  the  diagonal  element  of  Bardeen’s  transition  operator, 

Tkk  =  ( a  *  ek  Gvv* 


and  is  the  energy  level  calculated  as  if  the  barrier  were  very  wide.  With 
appropriate  choices  of  phase  for  the  (x),  O^(x),  one  naturally  defines 

The  0^(x)  are  large  only  in  the  right-hand  region,  but  there  is  a  long  tail 
(proportional  to  T^)  in  the  "wrong”  region.  Defining  the  a  and  b  operators 
according  to  (2),  the  Hamiltonian  becomes 


-  Z  ;+6*  +  VV)  + 


(5) 


We  have  not  bothered  to  write  out  the  contributions  from  the  interaction  term. 
This  will  be  taken  up  below. 

Neglecting  the  interaction  terms  then^we  can  calculate  the  time  derivative 
of  the  number  of  particles  in  right-hand  states .  We  find 


Bjj  =  (Nl°  -  Nr°)  T  sin  2Tt 


(6) 
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We  have  assumed  that  T  is  slowly  varying  with  k,  and  has  average  value  T.  We 

KK 

have  also  supposed  that  is  constant,  and  that  at  t  =  0,  the  numbers  of 

particles  present  in  the  right  (left)  states  is  N  . 

R(L) 

Formula  (6)  exhibits  the  oscillatory  behavior  well-known  from  elementary 
quantum  mechanics.  Such  behavior  is  not  to  be  expected  in  the  context  of  real 
metals,  where  one  would  anticipate  finding  a  constant  transition  rate  for  not- 
too-long  times.  If  (6)  does  not  show  a  constant  transition  rate  at  any  time,  It 

can  only  be  a  result  of  the  two  clever  definition  of  the  states  &  .  We 

* 

have  put  electrons  into  the  long  tail  of  the  on  the  "wrong"  side  of  the 


barrier,  already  accomplishing  the  greater  part  of  that  tunnelling  which  usually 
proceeds  at  a  constant  rate.  If  we  look  at  the  number  of  electrons  in  right-hand 
states  we  will  not  obtain  a  true  picture  of  the  tunnelling,  at  least  not  without 
further  arguments.  By  the  same  token,  one  must  handle  gingerly  any  states  which 
are  not  confined  rather  strictly  to  one  side  or  the  other. 


■W* 


IV 

The  other  alternative  is  to  find  a  separation  of  the  Hamiltonian,  as  in 
Eq.  (1) ,  based  on  states  which  are  really  confined  to  one  side  or  the  other.  A 
natural  set  of  states  is  given  by  the  eigenstates  of  H^  ,  where 

H,  "  1/,1-X) 


(7) 


and 


(V)  r  WO  , 

'  K  ,*<-&* 


(8) 


This  particular  choice  is  convenient,  but  not  mandatory.  We  call  the  eigenstates 
(J^(x)  and  also  define 

~  (9) 

Each  of  the  sets  is  complete  in  the  entire  space  (a  >|x|  ).  However, 

to  represent  a  function  which  has  large  amplitude  on  the  left  by  ,  it  will  be 

necessary  to  employ  large  k  values.  Let  us  introduce  0^  and  el#  operators  accord¬ 
ing  to 


=  i  ck  <4«-> «  r  «4  x.ui 


(10) 


We  next  employ  the  representation 

-  r  (  <t>„ t>)  +  X*W) 


Then 


=  Z  ( <&,pr)  -  Xj*)) 
/**-  ^  ~  ?  <t>±y  0^)  T  ( 


(11) 

(12) 

(13) 

(14) 
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We  have  next  to  compute  (  (j)^  This  matrix  element  is  diagonal  in  all  quantum 
numbers  except  the  momentum  component  conjugate  to  x.  We  shall  suppress  these  other 
quantum  numbers,  as  we  have  up  to  now. 

The  largest  matrix  elements  between  states  of  practically  the  same  energy 
are  found  to  be,  after  a  tedious  calculation, 


(15) 


(16) 


(17) 


kk 


where  we  have  introduced 

The  dots  represent  terms  of  higher  order  in  T. 

In  addition,  there  are  large  matrix  elements  (  )  when  k  is  so  large 

that  the  wave-length  on  the  left  Is  practically  given  by  m.  In  that  region,  however, 
we  find  (  (V  >  —  l'1 

Thus  we  find  to  sufficient  approximation,  that 


--  U,1 

/W  r  %t/(l t-CO  +  i?*.) 

and  there  are  no  large  matrix  elements  to  high  energy  states. 


(18a) 


(18b) 
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In  the  remainder  of  this  section,  we  shall  neglect  all  inter-electron 
interaction.  Substitution  of  expressions  (10)  -  (12)  and  (18)  into  the 
Hamiltonian  yields  at  once 

H  =  Z 

'  +V* +  T^)  A- 

*  I  fk(£*<c t  +  444) 

'hi  ^  c£ )  -  [Cfc"^«  +  C^  Ct] 

^  ^  ^  (VM4C>^ 

tc.  »\  >1  ,  -  _ 


«  - 


The  last  term,  as  well  as  others  indicated  by  dots,  may  be  dropped  if  we  wish  to  keep  only 
the  terms  of  lowest  order  in  the  tunnelling  exponential.  The  desired  separation  of  the 
single-particle  Hamiltonian  has  thus  been  achieved  with 

^R.  ^  (20a) 

/4  s 


(20b) 
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However, 


Irs  (C*  ^S~  0 


(21) 


which  is  of  the  order  of  .  This  means  that  the  eigenstates  of  H  +  H  are 

R  L 

not  just  the  products  of  eigenstates  of  H^,  separately. 

Let  us,  for  example,  calculate  the  current  through  the  barrier^ J ,  where 


(22) 


It  may  readily  be  shown  that 


-  -tl  <V4' 

r  C  ku.' 

(2) 

by  the  methods  of  Bardeen.  We  propose  to  calculate 


(23) 


(24) 


J  ~  (  $(Oj  Top  $<-0 

where  ^  (t)  is  the  state  reducing  at  t  =  0  to  the  product  state 

§<«■>-  $e<iLio>  (25) 

being  an  eigenstate  of 

It  is  easy  to  verify  that  the  failure  of  H  to  commute  with  is 
just  compensated  by  the  term  containing £  in  (19).  Indeed 

(Hk+H J  §k§Jo>  --  §  H  §  10 ?  +  in  -l  1  c>  (J6> 

to  the  first  order  tunnelling.  In  consequence,  we  may  drop  that  term,  at  the  same 
time  assuming  f  Cr"^  Jgj  -O.  This  leads  to  the  usual  formula  for  the  current 
proportional  to  1^4 1  ' 
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V. 

.We  turn  next  to  the  problem  of  interactions.  We  shall  treat  here  for 
simplicity  only  short  range  interactions  of  simple  type.  (The  range  should  be 
small  even  in  comparison  with  the  widths  of  the  barrier.) 

Perhaps  the  easiest  way  to  write  down  the  terms  of  the  interaction  Hamiltonian 
is  to  employ  the  representation 

iUxns  Z  +  I  X. . 

7  v 

(2?) 

We  have  introduced  the  notation 

Here  is  given  by  (17)  with  the  understanding  that  the  matrix  elements 

connecting  states  of  widely  different  energies  are  to  be  dropped. 

We  may  immediately  write  down  the  interaction  terms 

* "  [  3 1  U/;;S  sWc  4  &  wrVs 

+f + 

+?w'p% 


(29) 
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Here  ,  we  have  written 

tgi  fs-  **  *r 

=/i  ^-3)  fib)  &'•>>?,'*> 


(30b) 


lA/jy-s  -  //  w( *"j )  [  $rWvi»  tyj)  X*i»J  Ooe) 

I  ^  ^  A-  ^s“.  Aw-  Aw 

™p<yn  ~  J  J  W(x-y)[  $p\y)  X^>  f  (30d) 

*r 


Because  of  the  distinction  between 


^6  and 


we  see  that  the  usual  interaction 


terms,  Vr°; ,  are  slightly  modified,,  However,  this  is  of  no  interest.  On  the  other 
(1) 

hand  ,  W  is  a  term  leading  to  tunnelling.  It  will  be  small  of  order  T  and  could 

be  quite  important  numerically  for  strongly  interacting  systems.  The  remaining  terms 

(2) 

are  of  second  order  in  T  and  should,  for  the  sake  of  consistency,  be  dropped.  W  is 

of  some  interest,  however,  in  that  it  shows  that  a  direct  tunneling  of  two  particles 

(1) 

is  possible,  even  apart  from  the  higher  order  contributions  of  T  and  W 

(o) 

We  still  have  to  take  into  account  the  fact  that  the  large  interaction  terms,  W  , 

on  the  first  line  (28)  do  not  commute  with  one  another.  It  is  not  difficult  to  show 

(1) 

that  this  can  be  accounted  for,  at  least  in  first  orderly  a  change  in  W  .  The 
result  is  that  if  we  replace  the  second  line  of  (29)  by 

j-  I  |  Wf(.rJ  4  "VCV  C^Cf  ) 


+  (<kf^  ^v-cr  -i  ^/c^CrXs)  w 

then  we  may  regard  the  c1 s  and  d1 s  as  commuting  operators.  Expression  (31)  must  be 

,  %  (1)  ' 

thought  of  as  operating  to  the  right.  In  (31) ,  W  is  the  same  as  expression  (30b) 

Z  iXf 

with  the  tilde  removed  from  the  functions.  W  is  the  same  as  (30b)  with  both 


X  % 


functions  replaced  by 


functions  respectively. 


13 


In  summary,  we  have  succeeded  in  the  aim  of  splitting  the  Hamiltonian 
as  in  Eq.  (1) 0  The  parts  H^_  and  do  not  quite  commute,  but  we  have  shown,  at 
least  in  first  order,  that  this  can  be  completely  accounted  for  by  a  change  in 
the  transition  operator  T.  The  transition  term  is  the  same  as  written  down  by 
Bardeen  for  free  particles.  We  have  explicitly  found  the  correction  to  this 
in  the  presence  of  interactions.  This  correction  might  well  be  quantitatively 
significant. 


The  author  would  like  to  thank  Professors  Falk,  Ferrell,  and  GELick  for 
stimulating  conversations . 
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